Calculation of the Electron Self Energy for Low Nuclear Charge 



Ulrich D. Jentschura, 1,2 '0, Peter J. Mohr/fl, and Gerhard Soff 2, i 
1 National Institute of Standards and Technology, Gaithersburg, MD 20899-0001, USA 
2 Institutfur Theoretische Physik, TU Dresden, Mommsenstrafie 13, 01062 Dresden, Germany 



Abstract 

We present a nonperturbative numerical evaluation of the one-photon electron 
self energy for hydrogenlike ions with low nuclear charge numbers Z = 1 to 
5. Our calculation for the 15 state has a numerical uncertainty of 0.8 Hz for 
hydrogen and 13 Hz for singly-ionized helium. Resummation and convergence 
acceleration techniques that reduce the computer time by about three orders 
of magnitude were employed in the calculation. The numerical results are 
compared to results based on known terms in the expansion of the self energy 
in powers of Za. 
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Recently, there has been a dramatic increase in the accuracy of experiments that measure 
the transition frequencies in hydrogen and deuterium This progress is due in part 

to the use of frequency chains that bridge the range between optical frequencies and the 
microwave cesium time standard. The most accurately measured transition is the 1S-2S 
frequency in hydrogen; it has been measured with a relative uncertainty of 3.4 x 1CT 13 or 
840 Hz. With trapped hydrogen atoms, it should be feasible to observe the 1S-2S frequency 
with an experimental linewidth that approaches the 1.3 Hz natural width of the 25 level 
0J4]]. Indeed, it is likely that transitions in hydrogen will eventually be measured with an 
uncertainty below 1 Hz |5]U. 

In order for the anticipated improvement in experimental accuracy to provide better 
values of the fundamental constants or better tests of QED, there must be a corresponding 
improvement in the accuracy of the theory of the energy levels in hydrogen and deuterium, 
particularly in the radiative corrections that constitute the Lamb shift. As a step toward a 
substantial improvement of the theory, we have carried out a numerical calculation of the 
one-photon self energy of the IS" state in a Coulomb field for values of the nuclear charge 
Z — 1, 2, 3, 4, 5. This is the first complete calculation of the self energy at low Z and provides 
a result that contributes an uncertainty of about 0.8 Hz in hydrogen and deuterium. This 
is a decrease in uncertainty of more than three orders of magnitude over previous results. 

Among all radiative corrections, the largest by several orders of magnitude are the one- 
photon self energy and vacuum polarization corrections. Of these, the larger and historically 
most problematic is the self energy. Analytic calculations of the electron self energy at low 
nuclear charge Z have extended over 50 years. The expansion parameter in the analytic 
calculations is the strength of the external binding field Za. This expansion is semi-analytic 
[i.e., it is an expansion in powers of Za and \n(Za)~ 2 }. The leading term was calculated 
in |7|. It is of the order of a (Za) \n(Za)~ 2 in units of m e c 2 , where m e is the mass of the 
electron. In subsequent work |]7|-p5| higher-order coefficients were evaluated. 

The analytic results are relevant to low-Z systems. For high Z, the complete one- 
photon self energy has been calculated without expansion in Za by numerical methods 
fl26H37l. However, such numerical evaluations at low nuclear charge suffer from severe loss 
of numerical significance at intermediate stages of the calculation and slow convergence in 
the summation over angular momenta. As a consequence, the numerical calculations have 
been confined to higher Z. 

Despite these difficulties, the numerical calculations at higher Z could be used together 
with the power-series results to extrapolate to low Z with an assumed functional form in 
order to improve the accuracy of the self energy at low Z fl30fl ; up to the present, this approach 
has provided the most accurate theoretical prediction for the one-photon self energy of the 



IS state in hydrogen [38 



However, this method is not completely satisfactory. The extrapolation procedure gives 
a result with an uncertainty of 1.7 kHz, but employs a necessarily incomplete analytic 
approximation to the higher-order terms. It therefore contains a component of uncertainty 
that is difficult to reliably assess. Termination of the power series at the order of a (Za) 6 
leads to an error of 27 kHz. After the inclusion of a result recently obtained in p5f for the 
logarithmic term of order a (Za) 7 ln(Za)~ 2 the error is still 13 kHz. 

A detailed comparison between the analytic and numerical approaches has been inhibited 
by the lack of accurate numerical data for low nuclear charge. The one-photon problem is 
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especially well suited for such a comparison because five terms in the Za expansion have 
been checked in independent calculations. The known terms correspond to the coefficients 
An, Aio, A50) A32 and A31 listed below in Eq. (|]). 

The energy shift AZ^se due to the electron self energy is given by 



AE^=*^m e c 2 F(Za) 



it rr 



(1) 



where n is the principal quantum number. For a particular atomic state, the dimensionless 
function F depends only on one argument, the coupling Za. The semi-analytic expansion 
of F(Za) about Za = gives rise to the following terms: 



F(Za) = An H Za ) + A m + (Za) A 50 + (Za 
x [A 62 \n 2 (Za)- 2 + A31 \n(Za)- 2 + G SE (Za 



(2) 



where GsE(Za) represents the nonperturbative self-energy remainder function. The first 
index of the A coefficients gives the power of Za [including the (Za) A prefactor from Eq. 
(H)], the second corresponds to the power of the logarithm. For the IS ground state, which 
we investigate in this Letter, the terms An and Aio were obtained in [0-0 . The correction 
term A50 was found in Hl4"Hl6|l . The higher-order corrections A32 and A31 were evaluated 
and confirmed in [p^-21 



The results are 
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The Bethe logarithm lnfco has been evaluated, e.g., in as lnfco = 2.984 128 555 8(3). 

For our high-accuracy, numerical calculation of F(Za), we divide the calculation into 
a high- and a low-energy part (see Ref. Except for a further separation of the low- 

energy part into an infrared part and a middle-energy part, which is described in |^T| and 
not discussed further here, we use the same integration contour for the virtual photon energy 
and basic formulation as in 



The numerical evaluation of the radial Green function of the bound electron [see Eq. 

p. 296) 



(A. 16) in [28]] requires the calculation of the Whittaker function W KjfM (x 



sec 



over a very wide range of parameters n, /1 and arguments x. Because of numerical cancella- 
tions in subsequent steps of the calculation, the function W has to be evaluated to 1 part 
in 10 24 . In a problematic intermediate region, which is given approximately by the range 
15 < x < 250, we found that resummation techniques applied to the divergent asymptotic 
series of the function W provide a numerically stable and efficient evaluation scheme. These 
techniques follow ideas outlined in [031 and are described in detail in [BT . 



For the acceleration of the slowly convergent angular momentum sum in the high-energy 
part [see Eq. (4.3) in |]| I, we use the combined nonlinear-condensation transformation [04] 
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This transformation consists of two steps: First, we apply the van Wijngaarden condensation 



transformation |45fl to the original series to transform the slowly convergent monotone input 
series into an alternating series |46|. In the second step, the convergence of the alternating 
series is accelerated by the 5 transformation [see Eq. (3.14) in [44]] ■ The 5 transformation 



acts on the alternating series much more effectively than on the original input series. The 
highest angular momentum, characterized by the Dirac quantum number k, included in the 
present calculation is about 3 500 000. However, even in these extreme cases, evaluation of 
less than 1 000 terms of the original series is required. As a result, the computer time for 
the evaluation of the slowly convergent angular momentum expansion is reduced by roughly 
three orders of magnitude. The convergence acceleration techniques remove the principal 
numerical difficulties associated with the singularity of the relativistic propagators for nearly 
equal radial arguments. These singularities are present in all QED effects in bound systems, 
irrespective of the number of photons involved. It is expected that these techniques could 
lead to a similar decrease in computer time in the calculation of QED corrections involving 
more than one photon. 

In the present calculation, numerical results are obtained for the scaled self-energy func- 
tion F(Za) for the nuclear charges Z = 1,2,3,4,5 (see Table 1). The value of a used in 
the calculation is a^ = 1/137.036. This is close to the current value from the anomalous 
magnetic moment of the electron [47]. 



1/a = 137.035 999 58(52) . 

The numerical data points are plotted in Fig. [I], together with a graph of the function 
determined by the analytically known lower-order coefficients listed in Eq. @. 

In order to allow for a variation of the fine-structure constant, we repeated the calculation 
with two more values of a, which are 

l/a> = 137.035 999 5 and l/a< = 137.036 000 5 . 

On the assumption that the main dependence of F on Za is represented by the lower-order 
terms in ([3|), the change in F(Za) due to the variation in a is 



dF(Za) 
da 



5a 



8a 

-2 Au — + 

a 



Z A 50 + O(aln 2 a) 



5a 



for a given nuclear charge Z. Based on this analytic estimate, we expect a variation 

F(Za > ) - F(Za ) « F(Za ) - F(Za < ) « -9 x 10~ 9 

for the different values of a. This variation is in fact observed in our calculation. E.g., for 
the case Z = 2 we find 

F{2a < ) = 8.528 325 061(1), 
F(2a ) = 8.528 325 052(1) and 
F(2a > ) = 8.528 325 043(1). 

This constitutes an important stability check on the numerics and it confirms that the main 
dependence of F on its argument is indeed given by the lowest-order analytic coefficients 
An and A 50 . 
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In addition to the results for F(Zato), numerical results for the nonperturbative self- 
energy remainder function Gse{Zo:o) are also given in Table 1. The results for the re- 
mainder function are obtained from the numerical data for F(Zao) by direct subtraction of 
the analytically known terms corresponding to the coefficients An, A i0 , A 50 , A32 and Asi 
[see Eqs. (0,0)]- Note that because the dependence of F on Za is dominated by the sub- 
tracted lower-order terms, we have at the current level of accuracy Gse(-Z«<) = GsE^Zao) = 
Gse(^«>)- The numerical uncertainty of our calculaton is 0.8 x Z 4 Hz in frequency units. 

A sensitive comparison of numerical and analytic approaches to the self energy can be 
made by extrapolating the nonperturbative self-energy remainder function Gse(Z«) to the 
point Za = 0. It is expected that the function Gsv(Za) approaches a constant in the limit 
Za — > 0. This constant is referred to as Gse(0) = Aso- In the analytic approach, much 
attention has been devoted to the coefficient Aso PH-|2"3j|. The correction has proven to be 
difficult to evaluate, and analytic work on Aso has extended over three decades. A step-by- 
step comparison of the analytic calculations has not been feasible, because the approaches 
to the problem have differed widely. An additional difficulty is the isolation of terms which 
contribute in a given order in Za, i.e. the isolation of only those terms which contribute to 
Aso (and not to any higher-order coefficients). 

In order to address the question of the consistency of Aso with our numerical results, we 
perform an extrapolation of our data to the point Za = 0. The extrapolation procedure is 
adapted to the problem at hand. We fit Gse to an assumed functional form which corre- 
sponds to Aso? An and A 70 terms, with the coefficients to be determined by the fit. We find 
that our numerical data is consistent with the calculated value Aso = —30.924 15(1) f24|,f48|. 
It is difficult to assess the seventh-order logarithmic term An, because the extrapolated 
value for A 71 is very sensitive to possible eighth-order triple and double logarithmic terms, 
which are unknown. We obtain as an approximate result A 71 = 5.5(1.0), and we therefore 
cannot conclusively confirm the result |25| 



/139 \ 
A 71 = n -In 2] =4.65. 

Since our all-order numerical evaluation eliminates the uncertainty due to higher-order terms, 
we do not pursue this question any further. 

The numerical data points of the function Gse{Zo) are plotted in Fig. || together with 
the value Gse(0) = Aso = —30.92415(1). For a determination of the Lamb shift, the 
dependence of Gse on the reduced mass m r of the system has to be restored. In general, 
the coefficients in the analytic expansion (|2|) acquire a factor (m r /m e ) 3 , because of the 
scaling of the wave function. Terms associated with the anomalous magnetic moment are 
proportional to (m r /m e ) 2 ||49|| . The nonperturbative remainder function Gse is assumed to 
be approximately proportional to (m r /m e ) 3 , but this has not been proved rigorously. Work 



is currently in progress to address this question [[50 



We conclude with a brief summary of the results of this Letter, (i) We have obtained 
accurate numerical results for the self energy at low nuclear charge. Previously, severe 
numerical cancellations have been a problem for these evaluations, (ii) For a particular 
example, we have addressed the question of how well semi-analytic expansions represent 
all-order results at low nuclear charge. Our numerical data is consistent with the value 
A 60 = —30.92415(1) p4j , filfl . (iii) Numerical techniques |44| have been developed that 
reduce the computer time for the problem by about three orders of magnitude. 
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The calculation presented here is of importance for the interpretation of measurements 
in hydrogen, deuterium and singly-ionized helium and for the improvement of the Rydberg 
constant, because of recent and projected progress in accuracy. In the determination of the 
Rydberg constant, uncertainty due to the experimentally determined proton radius can be 
eliminated by comparing the frequencies of more than one transition ||. We have shown 
that an all-order calculation can provide the required accuracy if suitable numerical methods 
are used. 
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TABLES 



TABLE I. Scaled self-energy function and nonperturbartive self-energy remainder function for 
low-Z hydrogenlike systems. 



F(Za ) and G m {Za ) 



z 


F(Zao) 


Gse(Zoo) 


1 


10.316 793 650(1) 


-30.290 24(2) 


2 


8.528 325 052(1) 


-29.770 967(5) 


3 


7.504 503 422(1) 


-29.299 170(2) 


4 


6.792 824 081(1) 


-28.859 222(1) 


5 


6.251 627 078(1) 


-28.443 472 3(8) 
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FIGURES 



FIG. 1. The self-energy function F(Za). The points are the numerical results of this work, 
the curve is given by the analytically known terms that correspond to the coefficients listed in 
Eq. (|). 
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FIG. 2. Results for the scaled self-energy remainder function GsE(Za) at low Z. 
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